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ABSTRACT 

We combine a complete sample of 113 pointed observations taken with the Rossi X-ray Timing Ex- 
plorer between 1996-1999, monitoring observations taken with the Ryle telescope and the Green Bank 
Interferometer, and selected observations with the Very Large Array to study the radio and X-ray prop- 
erties of GRS 1915+105 when its X-ray emission is hard and steady. We establish that radio emission 
always accompanies the hard-steady state of GRS 1915+105, but that the radio flux density at 15.2 GHz 
and the X-ray flux between 2-200 keV are not correlated. Therefore we study the X-ray spectral and 
timing properties of GRS 1915+105 using three approaches: first, by describing in detail the properties 
of three characteristic observations, then by displaying the time evolution of the timing properties during 
periods of both faint and bright radio emission, and lastly by plotting the timing properties as a function 
of the the radio flux density. We find that as the radio emission becomes brighter and more optically 
thick, 1) the frequency of a ubiquitous 0.5-10 Hz QPO decreases, 2) the Fourier phase lags between hard 
(11.5-60 keV) and soft (2-4.3 keV) in the frequency range of 0.01-10 Hz change sign from negative to 
positive, 3) the coherence between hard and soft photons at low frequencies decreases, and 4) the relative 
amount of low frequency power in hard photons compared to soft photons decreases. We discuss how 
these results reflect upon basic models from the literature describing the accretion flow around black 
holes and the possible connection between Comptonizing electrons and compact radio jets. 

Subject headings: black hole physics — stars: individual (GRS 1915+105) — X-rays: stars 



1. INTRODUCTION 

The microquasar GRS 1915+105 is one of the most 
interesting galactic sources of radio, infrared, and X-ray 
emission. It was first discovered as a transient X-ray 
source by Granat in 1992 (Castro-Tirado et al. 1992), 
and later was observed to emit highly relativistic radio 
jets with intrinsic velocities greater than 0.9c (Mirabel & 
Rodriguez 1994; Fender et al. 1999b). High optical ex- 
tinction has prevented studies of the binary companion 
of GRS 1915+105, and there is no known orbital period 
or mass function for the system. However, the X-ray lu- 
minosity of GRS 1915+105 greatly exceeds the Edding- 
ton luminosity for a neutron star (Greiner, Morgan, and 
Remillard 1996), and the other galactic sources of highly 
relativistic (> 0.9c) jets, GRO J1655-40 (Orosz & Bai- 
lyn 1997; Shahbaz et al. 2000) and V4641 Sgr (Orosz et 
al. 2001), have mass functions greater than the maximum 
mass for a neutron star, 3 M Q . GRS 1915+105 is therefore 
thought to be a black hole binary system. The properties 
of microquasars in the galaxy are reviewed by Mirabel & 
Rodriguez (1999). 

Although GRS 1915+105 is often noted for exhibiting X- 



ray, radio, and infrared emission which is variable on time 
scales of seconds to minutes (Greiner et al. 1996; Pooley 
& Fender 1997; Mirabel et al. 1998; Eikenberry et al. 
1998), the source also exhibits steady emission which is 
reminiscent of black hole candidates such as Cyg X— 1 and 
GX 339—4. Both the steady and variable emission from 
GRS 1915+105 exhibit two (Muno, Morgan, & Remillard 
1999; Rao, Yadav, & Paul 2000) or three (Belloni et al. 
2000) basic modes of X-ray emission. In the two-state 
description, the hard state resembles the very high state 
of canonical black hole binaries (Morgan, Remillard, & 
Greiner 1997), as it exhibits X-ray flux with a power-law 
spectrum above 40 keV, prominent thermal emission from 
the accretion disk below ~ 5 keV, and a strong (up to 
15% RMS) 0.5-10 Hz quasi-periodic oscillation (QPO). 
Conversely, the soft state of GRS 1915+105 exhibits little 
X-ray emission above 40 keV and no 0.5-10 Hz QPO, al- 
though a stationary 67 Hz QPO and several low frequency 
(< 0.1 Hz) QPOs are often observed during this state. 

Long-term radio and X-ray monitoring of 
GRS 1915+105 (Rodriguez et al. 1995; Foster et al. 
1996; Harmon et al. 1997; Pooley & Fender 1997; Bandy- 
opadhyay et al. 1998) has revealed that periods of hard, 
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steady X-ray emission are often accompanied by bright, 
optically thick radio emission. This has been dubbed the 
"plateau state" by Foster et al. (1996; see also Fender 
2001). The radio plateau emission in GRS 1915+105 re- 
sembles steady, optically thick radio emission observed 
coincident with hard X-ray emission in the black hole 
candidates Cyg X— 1 (Hjellming, Gibson, & Owen 1975; 
Brocksopp et al. 1999), GX 339-4 (Corbel et al. 2000), 
and Cyg X-3 (Waltman et al. 1995; McCollough et al. 
1999). The connection between the X-ray and radio emis- 
sion in GRS 1915+105 is illustrated in Figure 1, where 
we plot the intensity as a function of time as observed 
in soft X-rays by the All-Sky Monitor (ASM) aboard the 
Rossi X-ray Timing Explorer (RXTE: first panel), and at 
radio wavelengths by the Ryle telescope (15.2 GHz; third 
panel). We also plot the the ASM hardness ratio HR2 
(5-12 keV / 3-5 keV; second panel), and the radio spectral 
index a taken from monitoring data at 2.25 and 8.3 GHz 
from the Green Bank Interferometer (i.e. S v oc v a \ only 
observations where S v at both wavelengths are greater 
than 20 mJy are plotted). Bright radio plateau emission 
and hard-steady X-ray emission is evident on several oc- 
casions (e.g. MJD 50730-50750), yet the radio emission 
is sometimes a factor of 10 fainter during other hard X- 
ray states (e.g. MJD 50450-50550). Observations with 
the Very Long Baseline Array during the plateau state 
reveal that the optically thick radio emission originates in 
a compact jet on order 10 AU across (Dhawan, Mirabel, 
& Rodriguez 2000). 

The radio plateau state is of further interest because it 
is often preceded and/or followed by optically thin radio 
flares in GRS 1915+105 (Fender et al. 1999b), as are 
similar hard-steady states with optically thick radio emis- 
sion in GX 339-4 (Corbel et al. 2000) and Cyg X-3 
(McCollough et al. 1999). The optically thin emis- 
sion in GRS 1915+105 originates from relativistic ejections 
of synchrotron-emitting blobs detached from the central 
source (Mirabel & Rodriguez 1994; Fender et al. 1999b; 
Dhawan et al. 2000). In Figure 1, optically thin radio 
flares precede and terminate the radio plateau emission 
which occurred between MJD 50730-50750, while weaker 
flares with shorter duration terminate the long period of 
weaker radio emission from MJD 50450-50550. These 
events are associated with flares in the ASM light curve, 
which suggests that the most powerful jets are associated 
with instabilities in the inner accretion disk. Moreover, 
there have been several occasions in which weaker radio 
and infrared flares recur at ^30 minute intervals, and are 
preceded by distinctive dip and flare patterns in X-rays 
(Pooley & Fender 1997; Eikenberry et al. 1998; Mirabel 
et al. 1998). 

During the plateau state of GRS 1915+105, the radio lu- 
minosity is low compared to the large optically thin flares, 
but the power in the radio jets may still represent a sig- 
nificant fraction of the accretion energy inferred from the 
X-ray emission (Falcke & Bicrmann 1999; Fender 2001). 
Since this emission originates from a compact jet within 
10 AU of the central source and is present only during the 
hard X-ray state, it seems very likely that the compact jet 
is influenced by the properties of the inner accretion flow, 
and vice versa. This has led Fender et al. (1999a) to sug- 
gest that the hard X-rays originate from inverse Compton 



scattering of seed photons at the base of the compact jet 
of relativistic electrons. This idea motivates further efforts 
to compare the properties of the radio emission to the X- 
ray spectral and timing information during the hard state 
of GRS 1915+105. 

The timing properties of GRS 1915+105 provide a 
promising means of probing the accretion flow near the 
black hole. In particular, the 0.5-10 Hz QPO observed 
during the hard state appears to provide a link between 
the optically thick accretion disk and the Comptonizing 
electrons, because (1) the frequency of the QPO increases 
as the thermal flux from the disk increases (Chen, Taam 
& Swank 1997; Trudolyubov, Churazov, & Gilfanov 1999; 
Markwardt, Swank, & Taam 1999), (2) the spectrum of the 
QPO amplitude is hard, indicating an origin in the power 
law component of the spectrum (Morgan et al. 1997), 
and (3) the QPO is only present when the power-law com- 
ponent of the X-ray spectrum is strong, while it is ab- 
sent when the power-law component is weak (Muno et 
al. 1999). Additional timing characteristics, such as the 
Fourier phase lags and the coherence function, may fur- 
ther constrain the relationship between hard and soft X- 
ray components. Reig et al. (2000) have found that the 
signs of the phase lags change from positive to negative 
as the frequency of the 0.5-10 Hz QPO increases, which 
has provided impetus for others (e.g. Nobili et al. 2000) 
to develop models describing the structure of the Comp- 
tonizing corona and how it is related to the optically thick 
accretion disk. 

In this paper, we use a complete sample of "hard- 
steady" observations of GRS 1915+105 (referred to as the 
X state by Belloni et al. 2000) to examine in detail the con- 
nection between the hard, steady X-ray emission and the 
steady radio emission. First, we establish that the hard- 
steady X-ray states of GRS 1915+105 always exhibit de- 
tectable radio emission. However, the intensities of the 2- 
200 keV X-ray and 15.2 GHz radio flux are not correlated, 
so we proceed to examine whether there are more subtle 
relationships between the spectrum and timing properties 
of the X-ray emission and the intensity and spectrum of 
the radio emission. We investigate these correlations in 
three different manners: first, by describing in detail the 
properties of three characteristic observations, then by dis- 
playing the time evolution of the timing properties during 
periods of both faint and bright radio emission, and lastly 
by examining the timing properties as a function of the 
radio flux density. 

2. OBSERVATIONS AND DATA ANALYSIS 
2.1. Ryle Monitoring Observations 

The Ryle Telescope at the Mullard Radio Astronomy 
Observatory, Cambridge, is primarily used for microwave- 
background measurements, but during gaps in that pro- 
gram it has been used to monitor variable galactic and 
extra-galactic sources. It operates at 15.2 GHz, and uses 
a single set of linearly-polarized feeds. In good observ- 
ing conditions, the RMS noise in a 1 minute integration 
is about 3.5 mJy. Since the measurements are unbiased, 
the noise level decreases as the square root of the observ- 
ing time. Details of the observing and analysis routine are 
given in Pooley & Fender (1997). 

We have averaged the radio measurements in two ways 
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for the purposes of this paper. We have used the aver- 
aged daily flux density measured with the Ryle Telescope 
when displaying monitoring light curves of GRS 1915+105 
(e.g. Figure 1). When comparing the Ryle measurements 
to the properties of pointed RXTE observations, we have 
averaged data taken within 0.5 days before and after the 
RXTE observation in order to sample the radio properties 
which most closely correspond to the X-ray observation 
(see Table 1). 

2.2. Green Bank Interferometer Monitoring Observations 

Monitoring with the Green Bank Interferometer has 
been performed simultaneously at 2.25 and 8.3 GHz on 
a 2.4 km baseline with 35 MHz bandwidth and dual circu- 
lar polarization. Most scans are 10-15 minute integrations, 
and the calibration procedure is described by Foster et al. 
(1996). Random (one sigma) errors in the GBI data are 
flux dependent: 4 mJy (2 GHz) or 6 mJy (8 GHz) for 
fluxes < 100 mJy, 15 mJy (2 GHz) or 50 mJy (8 GHz) 
for fluxes ~ 1 Jy. In addition, systematic errors arc in- 
troduced by atmospheric and hardware effects, which we 
estimate may approach 10% at 2 GHz and 20% at 8 GHz, 
occasionally higher for extreme local hour angles. Flux 
density measurements below 20 mJy approach the noise 
level of the instruments, and therefore we do not include 
measurements below this level in our analysis. 

As with the Ryle data, daily GBI measurements have 
been averaged when producing monitoring light curves, 
and GBI data taken within 0.5 days of a pointed RXTE 
observation have been averaged when comparing individ- 
ual observations (Table 1). In both cases, we also have 
computed the spectral index a from the average flux in 
the 2.25 and 8.3 GHz bands. 

2.3. Very Large Array Observations 

The VLA is a multi-frequency, multi-configuration aper- 
ture synthesis imaging instrument, consisting of 27 anten- 
nas of 25 m diameter. The receivers at 1.42, 5.0, 8.45, and 
15.2 GHz, have been used on various occasions, with 2 ad- 
jacent bands of 50 MHz nominal bandwidth processed in 
continuum mode. The corresponding 1-sigma sensitivities 
in 10 minutes are 0.06, 0.06, 0.05, and 0.2 mJy respectively. 
Switching between multiple frequencies takes about 30 s. 
The array configuration is varied every 4 months to cycle 
between 4 sets, A, B, C, and D, with maximum baselines 
of about 36, 11, 3.4, and 1 km. 

Observations of GRS 1915+105 while in its radio faint 
state have been obtained from the public archive in or- 
der to constrain the radio spectrum. The observations we 
have found are listed in Table 2, along with mean fluxes 
from Ryle and GBI observations taken within a half day 
of the VLA observations. For all of the observations re- 
ported here, GRS 1915+105 was unresolved by the syn- 
thesized beam (see Table 2) in any configuration at any 
wavelength. Other sources in the field of view (mainly at 
1.4 GHz) are well separated in the images, so there is no 
confusion with GRS 1915+105. 

Calibration and imaging were carried out with stan- 
dard tasks in the NRAO AIPS (Astronomical Image 
Processing System) package. For all of these obser- 
vations, the phase calibrator was 1925+211. The pri- 
mary flux density calibrator was 3C286 (1328+307) for 



all observations except 1996 November 14, when 3C48 
(0134+329) was used. In practice, the flux density er- 
rors are not set by the RMS receiver (thermal) noise 
stated above, but by errors in the flux density scale, es- 
timated to be 3-5% of the measurement, and/or source 
variability, depending on the occasion. More details of 
the VLA are given in the Observational Status Summary 
(http://www.aoc.nrao.edu/vla/obstatus/vlas/vlas.html), 
and the data reduction is treated in formal detail in Tay- 
lor, Carilli, & Perley (1999) and in the VLA cookbook 
(http: / / www.cv.nrao.edu/aips/cook.html) . 

2.4. Pointed RXTE Observations 

The Rossi X-ray Timing Explorer (RXTE) consists of 
three instruments: the All-Sky Monitor (ASM; Levine et 
al. 1996), the Proportional Counter Array (PCA; Jahoda 
et al. 1996), and the High Energy Timing Experiment 
(HEXTE; Rothschild et al. 1998). Figure 1 shows the 
ASM light curve of GRS 1915+105 along with daily radio 
measurements from Ryle and the spectral index from the 
GBI. There is no apparent correlation between the radio 
emission and the ASM count rate, although there is some 
correlation with the ASM hardness ratio HR2 (5-12 keV 
/ 3-5 keV). In addition, there has been a long campaign 
of weekly observations with with the PCA and HEXTE, 
which we use to study the detailed timing and spectral 
properties of the hard emission from GRS 1915+105. We 
define intervals as "hard-steady" if the ratio of the mean 
PCA count rates in the 12-60 keV to the 2-12 keV bands 
exceeds 0.05 and the mean variability in the count rate (2- 
60 keV) at Is time intervals is less than 12% for an RXTE 
orbit (typically 3000 s). In Figure 1 we have marked the 
pointed RXTE observations of the hard-steady state with 
vertical bars above the ASM light curve, and we have listed 
them in Table 1. The small dots represent the remainder 
of the pointed observations. We have analyzed all 113 ob- 
servations of GRS 1915+105 in the hard steady state taken 
with RXTE between 1996 April 15 (MJD 50188) through 
1999 March 22 (MJD 51259), during the time interval for 
the PCA instrument gain setting known as "epoch 3" . In 
order to improve the signal-to-noise of the spectral and 
timing measurements whenever possible, if the variability 
in the mean count rate in 16 s bins is less than 5% during 
a single day, we have analyzed all of the data for that day 
together. Otherwise, we have grouped the data while re- 
maining within the 5% variability limit (16 s bins), with 
the smallest unit being a single RXTE orbit. 

For each time unit, we have integrated 128-channel en- 
ergy spectra ( "Standard2" ) data from PCUs and 1 of 
the PCA, and 64 channel (archive mode) spectra from 
both clusters of HEXTE. In order to calculate the flux 
from GRS 1915+105, we have fit the spectra with the 
phenomenological model described in Muno et al. (1999). 
The model consists of the sum of a multi-color blackbody, 
a power law (with a high-energy cutoff when necessary), 
and a Gaussian line (representing iron emission between 
5-7 keV) . Interstellar absorption is taken into account as- 
suming a column density equivalent to 6 x 10 22 cm 2 . A 
normalization factor is allowed to account for differences 
in the effective area of the PCUs and the HEXTE clus- 
ters which are not taken into account by the response 
matrices. The standard FTOOLS 4.2 background sub- 
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traction algorithm for bright sources has been applied to 
the data, and a systematic uncertainty of 1% has been 
added to each spectral bin. We have computed an X-ray 
flux for each observation, by summing the bolometric flux 
implied by multi-temperature disk model (Fbb = 1-08 x 
10 _11 A^ ap pCrTg 1 pp ), the flux in the power law between 2-200 

keV (J™° 1.60 x 10- 9 N r E- r+1 cxp[-~E/E cut ]dE), and the 

flux in the iron line (N Giluss / (V2TrW Gauss ) J °° cxp[0.5(F- 

E Gauss ) 2 /W* auss }dE). 

We also have created power density spectra (PDS) and 
cross spectra (CS) for each 256 s interval of data with 2~ 7 
s time resolution and four energy channels, using combina- 
tions of the binned and event modes for each observation. 
The four energy channels used are as follows: 2-4.3 keV 
(henceforth referred to as the "low energy" band), 4.3-7.8 
keV, 7.8-11.5 keV, and 11.5-60 keV (referred to as the 
"high energy" band). The PDS and CS have been aver- 
aged for the entire time unit (an RXTE orbit or an obser- 
vation, as defined above) and logarithmically re-binned. 
The PDS have been corrected for dead-time effects and 
Poisson noise (Morgan et al. 1997), and have been nor- 
malized to the fractional RMS squared per Hz. To quan- 
tify the features of the PDS, we have fit PDS from each 
energy band with a model consisting of a power law with 
two breaks and Lorentzians for any QPOs present in the 
PDS. We have searched for the 67 Hz QPO during each 
of these hard-steady observations using data with 2~ 13 s 
time resolution in a single energy channel, but we did not 
detect it. Typical upper limits range from 0.1-0.5% (1 a). 

The coherence function and phase lags have been cal- 
culated from the CS in the manner described in Vaughan 
& Nowak (1997). We have used the standard convention 
that "hard" or "positive" lags indicate that higher energy 
photons lag behind lower energy photons. The uncertain- 
ties have been calculated for the case of high signal power 
and high measured coherence, using an estimate of the 
dead-time-corrected Poisson noise (Morgan et al. 1997). 
If the signal power is smaller than the estimated noise, or 
if the coherence is less than the noise divided by the sig- 
nal, we do not plot either the phase lags or the coherence 
functions. 

Guided by the long-term monitoring light curves in Fig- 
ure 1, we have defined subsets of the hard-steady observa- 
tions based upon the radio emission from GRS 1915+105 
in order to facilitate the comparison between the ra- 
dio and X-ray emission. We have defined "radio 
plateau" hard-steady conditions for observations for which 
GRS 1915+105 is brighter than 20 mJy at 15.2 GHz, and 
exhibits a flat radio spectrum with index (measured from 
the daily average of the flux) a > —0.2. We also have 
defined "radio steep" conditions for observations when 
GRS 1915+105 is brighter than 20 mJy at 15.2 GHz, 
and exhibits a optically thin radio spectrum with index 
a < —0.2. Figure 1 demonstrates that the radio steep 
conditions typically represent the transition into and out 
of periods of radio plateau emission, and occur during the 
decays of large optically thin radio flares (see also Fender 
et al. 1999b; Dhawan et al. 2000). Our discussion 
of radio steep observations will be limited, as the opti- 
cally thin emission probably results from radio ejecta that 
could evolve independently of the instantaneous X-ray con- 
ditions on 0.5 day time scales. 



Finally, we have defined "radio-faint" conditions for 
those hard-steady X-ray observations for which the daily 
average of the radio flux at 15.2 GHz is less than 20 mJy, 
because for these observations the radio flux at the GBI 
frequencies is too low to reliably measure a spectral in- 
dex. The radio faint conditions occur with two variants — 
those observations for which the X-ray emission is bright, 
and those for which it is faint. The behavior of the timing 
properties with respect to the radio emission is similar for 
both variants, although the shape of the X-ray spectrum 
and the values of the timing properties (e.g. the QPO fre- 
quency and the phase lags) are distinctive, and warrant a 
separate discussion. 

3. RESULTS 

We first investigate 1) whether all of the hard-steady ob- 
servations are coincident with radio emission, 2) whether 
the X-ray and radio flux are correlated, and 3) how the 
radio conditions (plateau, steep, and faint) are distributed 
in a plot of radio vs. X-ray flux. Figure 2 displays a plot 
of the radio flux density at 15.2 GHz (log scale) versus 
the 2-200 keV X-ray flux for each of the hard-steady ob- 
servations. It is clear that all hard steady observations 
have radio flux densities greater than 2 mJy. Panel 3 of 
Figure 1 shows that the radio flux from GRS 1915+105 
does drop below 1 mJy (e.g. MJD 50825-50875), but such 
intervals are associated with the soft X-ray state. Thus, 
radio emission always accompanies the hard-steady X-ray 
state, although a note of caution is appropriate for those 
observations associated with variable, optically thin radio 
emission, as noted above. However, Figure 2 also shows 
the radio and X-ray fluxes are not correlated in the hard- 
steady state. Likewise, we have found that the total 50-100 
keV flux, the power law flux, and the thermal (multi-color 
disk) flux are not correlated with the strength of the radio 
emission (not shown). 

The radio conditions are represented in Figure 2 as fol- 
lows: radio plateau observations are plotted as triangles, 
radio steep as circles, and radio faint as x's. Observa- 
tions for which the GBI spectral index has not been mea- 
sured are plotted only with error bars. The plot can 
be roughly split into four quadrants. The radio plateau 
observations are clearly localized in the upper left quad- 
rant, with 5i5.2 > 20 mJy and F x < 4 x 10~ 8 erg cm~ 2 
s _1 . Some radio steep observations are also present in this 
quadrant. The upper right quadrant contains radio steep 
observations and observations for which no spectral index 
is available from the GBI. These latter observations were 
taken between MJD 50313-50320 (before the current GBI 
monitoring program), a time interval which is included 
in the analysis of Bandyopadhyay et al. (1998). Dur- 
ing this time span a large radio flare occurred terminating 
a plateau state, and the radio flux density at 1.4 GHz 
was consistently larger than that at 3.3 GHz, which sug- 
gests that these points should be classified as radio steep. 
Therefore, the upper right quadrant (S15.2 > 20 mJy and 
Fx > 4 x 10 -8 erg cm~ 2 s _1 ) appears to include only ra- 
dio steep observations, during which the radio and X-ray 
emission may vary independently since the radio measure- 
ments may be dominated by the evolution of an expanding 
jet detached from the central source. For this reason, we 
choose to ignore the upper right quadrant. 
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The lower quadrants, with 515.2 < 20 mJy, contain radio 
faint observations by construction. After examining the 
temporal clustering of these radio- faint observations, we 
feel that it is useful also to consider them in two quadrants: 
an X-ray faint quadrant with Fx < 4 x 10~ 8 erg cm~ 2 
s" 1 , and an X-ray bright quadrant with Fx > 4 x 10~ 8 erg 
cm~ 2 s _1 . This separation is further motivated by their 
differing energy spectral and timing properties, which are 
described in the next section. 

3.1. Characteristic Observations 

From each of the three quadrants under consideration 
from the radio versus X-ray flux plane (Figure 2), we have 
selected a characteristic observation in order to illustrate 
the range of X-ray spectral and timing properties within 
the hard-steady X-ray state. In Figures 3, 4, and 5, the left 
panel is a radio faint, X-ray faint observation (MJD 50488; 
1997 February 9); the center panel is a radio plateau, X- 
ray faint observation (MJD 50737; 1997 October 16); and 
the right panel is a radio faint, X-ray bright observation 
(MJD 50708; 1997 September 17). The radio and X-ray 
fluxes from these three observations are indicated in Fig- 
ure 2. 

3.1.1. X-ray Energy Spectra 

Figure 3 illustrates the results of the spectral fits for the 
three representative observations with differing values of 
the radio and X-ray flux. The spectrum in the left panel 
is from a radio faint observation when the X-ray emis- 
sion was also relatively faint (2.2 x 10~ 8 crg cm~ 2 s -1 ). 
The soft emission during this observation can be modeled 
(xt ~ 1) using a cool multi-temperature blackbody with 
an apparent inner disk temperature of T app ~ 0.8 keV and 
an apparent radius of i? app ~40-80 km (assuming a dis- 
tance of 11 kpc and an inclination angle of 66°; see Fender 
et al. 1999b). The power law has a photon index r ~ 2 
and requires an exponential cutoff at about 70 keV. The 
results of this spectral fit are characteristic of observations 
during which both the X-ray and radio flux were low. 

The spectrum in the middle panel is from an observation 
which exhibited radio plateau conditions when the X-ray 
emission was also faint (2.6 x 10~ 8 erg cm~ 2 s _1 ). Despite 
the similarity in X-ray flux to the left panel, the shape of 
the X-ray spectrum is quite different. The thermal compo- 
nent of the model appears hot (T app > 3 keV) with a small 
apparent inner radius (-R app < 5 km). This quasi-thermal 
component is not consistent with simple models of reflec- 
tion from an accretion disk (e.g. Magdziarz & Zdziarski 
1995). It is possible that this broad thermal bump rep- 
resents emission from a relatively small area of the disk 
which has been heated by, for example, a magnetic flare 
(di Matteo, Celotti, & Fabian 1999) or a spiral shock wave 
(Tagger & Pellat 1999), or that the disk develops an op- 
tically thick atmosphere which Comptonizes the emergent 
spectrum, so that the spectral hardening and decreased 
apparent radius are illusory effects of increased scattering 
in a "puffed up" disk (Merloni, Fabian, & Ross 2000). The 
photon index of the power law is T ~ 2.7, and a cut-off at 
about 80 keV is required for observations with longer inte- 
gration times and hence better statistics at energies above 
100 keV. We note that Rao et al. (2000a) have presented 
a different model for spectra of GRS 1915+105 such as 



these, consisting of a multi-temperature disk, a thermal- 
Compton spectrum (Sunyaev & Titarchuk 1980), and a 
power law. For the purposes of this paper, we prefer to 
use our simpler model with fewer free parameters. 

The right panel is from an observation when the radio 
emission was faint (6 mJy) and the X-ray emission was 
bright (8.4 x 10~ 8 erg cm~ 2 s _1 ). The spectrum is similar 
to the radio bright, X-ray faint observation in the center 
panel, except that the relative contribution of the power 
law flux increases dramatically. The thermal component 
is hot (Tapp > 2.5 keV) and has a small apparent inner 
radius (i? app < 10 km), while the power law index is steep 
(r ~ 3.2). The shape of the energy spectrum on 1997 
September 17 is characteristic of all X-ray spectra from 
the hard state of GRS 1915+105 when the total flux is 
greater than ~ 4 x 10~ 8 erg cm~ 2 s _1 . The count rate 
from GRS 1915+105 tends to vary on time scales of hours 
when the source is bright in X-rays, so we only integrate 
spectra for individual RXTE orbits. With these exposure 
times, we can measure the power law only up to about 
100 keV. However, when we combine spectra with similar 
photon indices, we find that the data are consistent with 
no cut off in the power law below 200 keV. 

The changes in the X-ray spectrum clearly do not track 
the radio flux from GRS 1915+105 in a monotonic man- 
ner. The power law component steepens dramatically as 
the radio emission increases in the X-ray faint state (left 
and center panels of Figure 3), but is steepest when the X- 
ray emission is bright and the radio emission is faint (right 
panel). At low X-ray fluxes (< 4 x 10~ 8 erg cm~ 2 s _1 ), 
the thermal component is hot and has a small apparent 
inner radius during radio plateau observations, while it is 
cool and has a large apparent inner radius when the radio 
emission is faint (see also Muno et al. 1999). 

3.1.2. Power Density Spectra and Cross Spectra 

The power spectra (here displayed as the frequency 
times the power density), Fourier phase lags, and coher- 
ence functions of the three observations from Figure 3 are 
displayed in Figure 4. In the top panel, two lines are plot- 
ted for each power spectrum: the black line indicates the 
2-4.3 keV (soft) energy band, while the the grey line indi- 
cates the 11.5-60 keV (hard) band. The PDS in all three 
examples of Figure 4 have relatively flat spectra with a 
low frequency break, a QPO between 0.5-10 Hz, and a 
dramatic break at frequencies higher than the QPO. In 
the middle panel of Figure 4, we plot the phase lags be- 
tween the 11.5-60 keV and the 2-4.3 keV bands, and in 
the bottom panel we plot the coherence function for the 
same two energy bands. To simplify the comparison be- 
tween the power spectra, the phase lags, and the coher- 
ence functions, we have indicated the frequencies of the 
QPOs in our fits to the 2-4.3 keV band with vertical black 
short-dashed lines, and the continuum break frequencies 
with vertical black long-dashed lines. We have also indi- 
cated the break frequency in the 11.5-60 keV band with 
grey long-dashed lines (see below). For the following dis- 
cussion we define three frequency regions related for the 
features of the power spectrum: region 1 represents fre- 
quencies lower than the low frequency continuum break, 
region 2 frequencies between the low frequency break and 
the QPO frequency (less its width), and region 3 frequen- 
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cies higher than the QPO frequency (or its harmonic if 
present). 

The first point of interest is that the QPO fre- 
quency is not strictly correlated with the X-ray flux from 
GRS 1915+105. From the spectra in Figure 3, it is clear 
that the X-ray flux increases from left to right, while in 
Figure 4 the QPO frequency is lowest in the center panel 
(radio plateau conditions). In the following sections, we 
will demonstrate that the QPO frequency is always lowest 
during radio plateau conditions. 

The second and most obvious difference between the 
timing properties of these three observations is in the 
Fourier phase lags. In region 2 (between the low frequency 
break and the QPO), the phase lags are negative when 
the radio emission is faint (left and right panels) , and pos- 
itive during radio plateau emission (center panels; compare 
Reig et al. 2000). This behavior is in sharp contrast to 
sources such as GX 339—4 and Cyg X— 1, for which the 
phase lags are always positive. Only the phase lags dur- 
ing the radio plateau state resemble those of GX 339—4 
(Nowak et al. 1999c) and Cyg X-l (Cui et al. 1997). 

The third notable result concerns the coherence function 
in region 1 (i.e. frequencies lower than the first break). 
The coherence function between the 11.5-60 keV and the 
2-4.3 keV bands is lowest (~ 0) during the radio plateau 
conditions (center panel), somewhat higher (~ 0.4) when 
both the radio and X-ray emission are faint (left panel), 
and nearly unity when the radio emission is faint and the 
X-ray emission is bright (right panel). Only the latter 
case is reminiscent of Cyg X— 1 in its hard state (Cui et 
al. 1997; Nowak et al. 1999a). 

Finally, a comparison of the amplitude and strength of 
the power spectra as a function of energy reveals further 
differences in each of the three observations. While the ra- 
dio emission is faint, (right and left panels) the root mean 
power at low frequencies is greater in the 11.5-60 keV en- 
ergy band (gray curve) than in the 2-4.3 keV band (black 
curve). However, during radio plateau conditions (center 
panel) the root mean power is nearly equal in all energy 
bands. The shape of the power spectrum also changes — 
the low frequency break is strongest in the low energy band 
in the radio plateau observation (center panel), while the 
break is most apparent in the high energy band in the ra- 
dio faint, X-ray faint observation (left panel). When the 
X-ray emission is bright, the low frequency break is evident 
in both energy bands (right panel). 

There are several other interesting details in the timing 
properties presented in Figure 4, although they are not 
directly relevant to the question of how the timing prop- 
erties are related to the radio emission. The harmonic of 
the QPO tends to be strongest (~ 10%) when count rate is 
lowest (left and center panels), but is not apparent during 
X-ray bright observations (right panel). A weak feature 
at approximately half of frequency of the strongest QPO 
peak is also evident during observations when both the 
X-ray and radio flux are faint (right panel), with an RMS 
amplitude of about ~1% and Q values of 1-10. 

Finally, a weak QPO is evident between 0.01-0.04 
Hz during a stretch of radio faint observations from 
MJD 50450-40560 (left panel), with an RMS amplitude 
of about 0.1% in the 2-4.3 keV energy band. This low 
frequency QPO disappears at higher energies. Its ccn- 



troid frequency decreases as the frequency of the 0.5-10 
Hz QPO decreases (not shown). Since the 0.01-0.04 Hz 
QPO is present only in the low energy band, the coher- 
ence function at the QPO frequency is near zero. 

3.1.3. Relation of the Timing Properties to the Energy 

Spectrum 

The spectral fits which we have made to GRS 1915+105 
indicate that two spectral components are present (Fig- 
ure 3). Most models of the emission from X-ray binaries 
invoke inverse-Compton scattering to produce hard pho- 
tons from soft seed photons, and consequently the vari- 
ability in the hard and soft energy bands should be inti- 
mately related. However, before finalizing any conclusions 
regarding the phase lag and coherence measures and how 
these change with photon energy, we need to examine how 
these timing properties depend upon the amount of the 
flux which is contributed by each of the spectral compo- 
nents in the relevant energy bands. 

In Figure 5, we examine the dependence of the timing 
properties on photon energy for the same three represen- 
tative observations as in Figures 3 and 4 (diamonds con- 
nected by solid lines). In all of the panels, the fraction of 
the count rate contributed by the disk component is dis- 
played with triangles connected with dashed lines, using 
the scale displayed on the right axis. Here it can be seen 
that the thermal component contributes no more than half 
of the flux in all energy bands for the observations which 
we analyze in this paper. Examining next the top panels, 
we find that the presence of the thermal component does 
not greatly affect the RMS power in a given energy band. 
In the second panels, the slopes of the phase lags show no 
obvious breaks; we might expect changes in the slopes if 
the soft and hard spectral components separately affected 
the phase lags. Finally, while the fraction of the flux in 
the thermal component increases in a similar manner dur- 
ing the X-ray faint, radio plateau (center panel) and the 
X-ray faint, radio bright observation (left panel), the co- 
herence function at low frequencies plummets dramatically 
with increasing energy in the plateau observation (center 
panel), but is near unity when the X-ray flux is bright 
(right panel) . We conclude that the changes in the timing 
properties arc not related to the fraction of the emission 
which can be modeled with a thermal component. 

3.2. Time Evolution of X-ray and Radio Properties 

The evolution of the parameters of the power spec- 
trum, phase lags, and coherence function can be viewed 
as a function of time during the long periods in which 
GRS 1915+105 remains in the hard-steady state. The re- 
sults demonstrate that the three observations selected for 
Figures 3-5 are representative of the groups which we have 
defined. In Figure 6 we show two time intervals contain- 
ing the hard-steady states of GRS 1915+105. In the first 
the radio emission is faint (MJD 50400-50600, left pan- 
els), and the second coincides with strong radio plateau 
emission (MJD 50725-50755, right panels). In the top 
three panels from both time intervals we plot: a) the in- 
tensity of GRS 1915+105 as viewed by the ASM, b) the 
radio flux measured with the Ryle telescope at 15.2 GHz, 
and c) the radio spectral index between 2.25 and 8.3 GHz 
(a) measured with the GBI. Notice that before the period 
of radio-plateau emission (right panels) the radio spectral 
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index increases smoothly and becomes positive (more op- 
tically thick), and then it decreases symmetrically after- 
wards. Optically thin radio flares that are much more 
prominent at 2.25 and 8.3 GHz (not shown) occur both 
before and after this interval. 

In panel d), we plot the frequency of the 0.5-10 Hz QPO. 
Trudolyubov et al. (1999) and Chen et al. (1997) sepa- 
rately established that the QPO frequency tracks the flux 
and the count rate in the PCA during hard-steady states 
similar to those displayed. However, it is clear that QPO 
frequency is systematically higher during radio faint emis- 
sion (> 2 Hz) than during radio plateau emission (< 2 Hz), 
despite the fact that the ASM count rate is lower during 
the faint radio emission. 

In panel e) of Figure 6 we plot the average phase lags 
in region 2 between the 2-4.3 keV band and 1) the 7.8- 
11.5 keV band (diamonds connected by the solid line), 
and 2) the 11.5-60 keV band (triangles connected by the 
dashed line). During radio-faint hard-steady conditions, 
the phase lags are consistently negative and have a value 
of about —0.2 radians. During radio plateau emission, the 
phase lags are positive and increase strongly with energy. 
Only during radio plateau conditions are the X-ray phase 
lags positive, as might be expected from a simple model 
which seeks to explain the timing properties using Comp- 
tonization (see Section 3). 

In panel f) of Figure 6 we plot the average coherence 
function in region 1 between the 2-4.3 keV band and 1) 
the 7.8-11.5 keV bands (diamonds connected by the solid 
line), and 2) the 11.5-60 keV bands (triangles connected 
by the dashed line) as a function of time. The coherence 
over this frequency range drops dramatically when the ra- 
dio emission is strong and optically thick, as shown in the 
representative observation of the radio plateau state in the 
center panel of Figure 4. 

In the bottom panel we plot the total of the RMS noise 
at frequencies lower than the QPO (region 1 and region 2 
combined). The absolute values of the RMS low frequency 
noise should not be compared directly from observation to 
observation, because as the QPO frequency varies, so does 
the upper limit to the frequencies which we consider. How- 
ever, we have plotted the RMS low frequency noise for two 
energy bands: the 2-4.3 keV band (diamonds connected 
by a solid line) and the 11.5-60 keV band (triangles con- 
nected by dashed lines). It is clear that there is much 
more low frequency power in the 11.5-60 keV band than 
in the 2-4.3 keV band during radio-faint hard X-ray condi- 
tions, while during radio-plateau conditions, the amount 
of power in both energies is comparable. This can also 
be seen by comparing the radio plateau and X-ray- faint, 
radio-faint observations in Figure 4. 

3.3. Timing Properties and the Radio States 

We now examine directly the relationship between the 
radio and the X-ray properties of GRS 1915+105. In Fig- 
ures 7 and 8, the various symbols denote the emission 
states in Figure 2. 

In Figure 7 we plot as a function of the radio flux den- 
sity at 15.2 GHz: a) the QPO frequency, b) the phase lag 
between the 2-4.3 keV and 11.5-60 keV bands in region 
2, c) the coherence for the same energy bands in region 
1, and d) the ratio of the low- frequency (regions 1 and 2) 



power in the 11.5-60 keV to 2-4.3 keV energy bands. At 
first glance, there appear to be continuous distributions in 
the radio and timing properties in Figure 7, but in fact the 
populations of radio faint observations and radio plateau 
observations are well-separated, with the radio steep ob- 
servations falling in between the two. Moreover, the x-axes 
of the panels are plotted with a logarithmic scale, which 
emphasizes the large range in the radio flux for the various 
hard-steady observations. 

Four conclusions can be drawn from Figure 7 regarding 
the changes in the timing properties as a function of the 
radio flux. The QPO frequency tends to decrease with in- 
creasing radio flux, and almost all of the observations with 
radio plateau conditions (triangles) have QPO frequencies 
less than 2 Hz. Second, during radio plateau hard-steady 
emission the hard photons lag the soft, while during ra- 
dio faint (x's) conditions the hard photons lead the soft. 
Third, the coherence at low frequencies is lowest during 
radio plateau emission. Finally, the ratio of low frequency 
power in 11.5-60 keV photons to the power in 2-4.3 keV 
photons is nearly unity only when the radio flux is largest, 
during radio plateau conditions. 

Using our large dataset which incorporates a variety of 
radio conditions, we next examine correlations which other 
authors have reported to be well-defined. Two studies 
(Chen et al. 1997; Trudolyubov et al. 1999) have ex- 
amined the frequency of the 0.5-10 Hz QPO as a function 
of the PCA count rate (MJD 50275-50333) and the 3-20 
keV X-ray flux (MJD 50363-50563), respectively. Both 
have found significant correlations among these values. In 
the top panels of Figure 8 we plot the QPO frequency as 
a function of total X-ray flux (defined in Section 2.3; left 
panel), thermal flux (middle panel), and power law flux 
(right panel). We confirm that the frequency of the 0.5- 
10 Hz QPO is tightly correlated with the X-ray flux from 
GRS 1915+105, but we also find that the correlations are 
not single- valued. In the left panel there exists a signif- 
icant population of points from radio faint observations 
that have QPO frequencies between 2-5 Hz, for which the 
X-ray fluxes are low compared to those expected from the 
general trend. 

A plot of the average phase lags between 2-4.3 and 11.5- 
60 keV photons in region 2 as a function of the total, ther- 
mal, and power law fluxes (bottom panels of Figure 8) 
reveals that there is no correlation between the X-ray flux 
and the sign or magnitude of the phase lags. The posi- 
tive phase lags therefore appear to occur independently of 
any given X-ray flux level, and are better predicted by the 
presence of optically thick radio emission. 

3.4. Characteristics of the Faint Radio Emission 

We can use the archival VLA observations to make some 
conclusions about the nature of the faint radio emission. 
The multi-frequency radio observations of GRS 1915+105 
are listed in Table 2. The first two observations occur 
within a day of radio faint, X-ray bright observations. In 
both cases, the VLA observed significant flaring, while 
Ryle observations, which were more nearly simultaneous 
with the RXTE observations, indicate that the radio emis- 
sion was near 5 mJy within 12 hours of the PCA obser- 
vations. It appears that the radio emission is unsteady 
within several hours of the X-ray bright hard-steady state, 
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and we can not rule out that the X-ray emission is also un- 
steady on these time scales. Since the VLA observations 
show clear evidence for flaring while the Ryle observations 
do not, we can not infer the radio spectrum of the emission 
during the X-ray bright, hard-steady state. 

The remainder of the observations in Table 2 are X-ray 
faint. There is clear evidence in the GBI data for ~30 
mJy optically thin flares on 1996 December 28, 1997 Jan- 
uary 12, and 1998 October 08. The flares decay on time 
scales of several hours, and no RXTE observations were 
taken within 0.5 day of these flares. Since GBI flux mea- 
surements below 20 mJy are uncertain, there remain two 
observations during which the VLA observations demon- 
strate that the radio emission was faint and steady. On 
1998 September 14, 15.2 GHz Ryle and 5.0 GHz VLA ob- 
servations separated by 6 hours provide a spectral index 
a = —0.5 ±0.1, which suggests that the emission is opti- 
cally thin. However, it is certainly possible that the radio 
flux varied by a few mJy during those several hours. Dual 
frequency VLA observations were made at 5.0 and 15.2 
GHz on 1998 September 29, which indicate a = 0.12±0.01. 
This suggests that the emission is optically thick, similar 
to a weak version of the radio plateau state. However, 
with only a single secure measurement of a during the ra- 
dio faint state, more observations are certainly warranted. 

4. DISCUSSION 

We begin the discussion with a summary of our three 
hard-steady conditions: 

• Radio Faint (e.g. MJD 50450-50550; MJD 50708) 
The radio emission has a mean flux density of about 
5 mJy at 15.2 GHz. When the X-ray flux is faint 
(< 4 x 10~ 8 erg cm~ 2 s _1 ), the X-ray spectrum is 
best described as the sum of a cool (T app <2 keV, 
R&pp ~ 60 km) multi-temperature blackbody and 
a power law of index L ~ 2.1 with a cut-off near 
70 keV. When the X-ray emission is brighter (> 
4 x 10~ 8 erg cm~ 2 s _1 ), the X-ray spectrum is best 
described with a hot (>2 keV) multi-temperature 
black body and a steep power law (T ~ 3.2) with- 
out a cut off below 200 keV (see Section 2.2). The 
QPO frequency is greater than 2 Hz, phase lags are 
negative (or, soft, indicating that hard photons pre- 
cede soft photons) at intermediate continuum fre- 
quencies, the coherence at low frequencies in the 
continuum is ~ 0.5-1.0 (11.5-60 keV compared to 
2-4.3 keV), and there is much more low frequency 
power in the high energy band than the low energy 
band. 

• Radio Plateau (e.g. MJD 50730-50748) The ra- 
dio emission is bright (~100 mJy at 15 GHz) and 
optically thick (a > —0.2). The X-ray spectrum 
is best described by a hot multi-temperature black 
body (T^pp > 3 keV and R app < 5 km) and a power 
law (r ~ 2.5) which could be cut-off at E > 80 keV. 
The QPO frequency reaches its lowest values (< 2 
Hz), phase lags are positive (hard photons lag soft 
photons), the coherence is low (particularly at low 
frequencies), and there is about equal broad band 
RMS power at low frequencies in all energy bands. 



• Radio Steep (e.g MJD 50727-50730) These obser- 
vations represent the transition into and out of ra- 
dio plateau conditions. The radio flux is bright, 
while the radio spectrum is optically thin (a < 
—0.2). The X-ray flux is often higher than during 
radio plateau conditions, but the spectrum is still 
best described by a hot multi-temperature black 
body (T^pp > 3 keV and i? app < 10 km) and a 
steep power law (L ~ 2.5) with a cut-off above 
80 keV. The timing properties tend to be similar 
to radio faint conditions (although there are excep- 
tions): the QPO frequency is between 2-5 Hz, the 
phase lags are negative, the coherence is around 0.5, 
and there is more broad band power at low frequen- 
cies in the highest energy band. 

Previous observations have established the nature of the 
bright radio emission from GRS 1915+105. The optically 
thick emission during plateau conditions has been resolved 
as a compact jet of relativistic electrons (Dhawan et al. 
2000), similar to ones which may be present in Cyg X— 1 
(Stirling, Spencer, & Garret 1998) and GX 339-4 (Cor- 
bel et al. 2000). Our observations of the faint radio emis- 
sion provide some evidence that it is similar to a weak 
radio plateau state. Radio flares with steep spectra origi- 
nate from material which has been ejected from the central 
source (Fender et al. 1999b; Dhawan et al. 2000), and 
so is most likely decoupled from the instantaneous con- 
ditions of the X-ray emitting regions. In the following 
discussion we will focus on comparing the radio faint and 
radio plateau conditions, during which the radio and X-ray 
properties of GRS 1915+105 are closely related. 

The timing properties of GRS 1915+105— the QPO 
frequency, phase lags, coherence function, and broad band 
noise — appear to be highly dependent upon the radio con- 
ditions (Figure 7). However, it is not immediately obvious 
whether the changes in the timing properties are part of 
the mechanism which produces a compact jet, or whether 
they result from changes in the accretion flow induced 
when the compact radio jet is strong. Since the location 
and origin of the Comptonizing electrons which generate 
hard X-rays is unknown, the possible connection between 
these and the jet is doubly uncertain. Nevertheless, we 
may explore interpretations of these results using simple 
models for the accretion flow under consideration in the 
literature. 

In Section 4.1, we consider the extent to which Comp- 
tonization of a soft input signal by a static corona at the 
base of a compact radio jet is consistent with the phase lags 
and X-ray spectrum of GRS 1915+105. In Sections 4.2- 
4.3, we examine two generic models which incorporate 
more elements of the accretion flow. In both models, cool 
(T <~ 1 keV) thermal emission originates from an optically 
thick accretion disk, while the hot power law component 
originates from inverse Compton scattering. The first ba- 
sic model places the Comptonizing electrons in a spherical 
corona within the inner radius of the accretion disk, while 
the second geometry assumes that the relativistic electrons 
are part of a planar corona sustained by magnetic flares 
above the disk. 

4.1. Generic Comptonization 
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Hard phase lags are generic properties of the low-hard X- 
ray state of black hole candidates such as Cyg X— 1 (Nowak 
et al. 1999a) and GX 339-4 (Nowak et al. 1999c). Sur- 
prisingly the phase lags in GRS 1915+105 change from 
negative to positive as the strength of the radio emission 
increases (Figure 7b). Since Fender et al. (1999a) have 
suggested that the Comptonizing corona in the hard state 
of black hole candidates represents the base of a compact 
radio jet, we feel it is important to explore the possibility 
that the negative phase lags during radio faint observa- 
tions of GRS 1915+105 are intrinsic to the region close to 
the black hole, while the hard phase lags generated during 
the radio plateau state are due to the Comptonization of 
the intrinsic input signal in a large corona at the base of 
the compact radio jet. 

The hard phase lags in GRS 1915+105 may constrain 
the structure of the Comptonizing region, because the 
phase is roughly constant as a function of frequency (Fig- 
ure 4). This is contrary to the constant time lags that 
would be expected from a homogeneous corona. Instead, 
Kazanas, Hua, & Titarchuk (1997) and Bottcher & Liang 
(1999) have demonstrated that constant phase lags can be 
produced in an isotropically illuminated, spherical Comp- 
ton cloud with a uniform temperature and a density which 
decreases with radius as since photons have an equal 
probability of scattering per decade of radius in such a 
cloud (see Nowak et al. 1999b for further discussion of 
these models). High-frequency variability originates from 
photons which have scattered from small radii, since the 
signal is washed out if it scatters with long time delays. 
On the other hand, the time lags of low-frequency signals 
are dominated by photons which have scattered over large 
radii. This produces time lags which decrease as a func- 
tion of frequency, or phase lags which are nearly constant 
as a function of frequency. 

Several authors have pointed out that the phase lags 
which are observed during the hard states of black hole 
candidates imply large scales for inverse Compton scat- 
tering, which raises the question of how energy is trans- 
ported to support a large, hot corona (Poutanen & Fabian 
1999; Nowak et al. 1999b). In the plateau state of 
GRS 1915+105 a phase lag of 0.5 radians at 0.3 Hz (mid- 
dle panel of Figure 4) would imply a light travel time of 
0.3 s, or a distance of Rt ~ 8 x 10 9 cm (where r is the 
optical depth of the corona). For a 10 M Q black hole, this 
represents thousands of Schwarschild radii. However, this 
distance is much smaller than the size of the compact jet 
observed by Dhawan et al. (2000), 10 AU = 5 x 10 14 cm, 
which is consistent with the assumption that the Comp- 
tonizing electrons reside at the base of the jet. If the 
Compton cloud is the base of a compact jet, this reduces 
the difficulty in explaining how energy is transported to 
support a very large corona, as the jet may carry up to 
10 39 erg s" 1 of power to much larger radii (Falcke & Bier- 
mann 1999). 

It is possible that the radio faint states are associated 
with a smaller Compton corona located near the inner ac- 
cretion disk. In this case, negative phase lags could be 
caused by any number of effects intrinsic to the accre- 
tion disk, such as waves similar to those of Nowak et al. 
(1999a) which propagate outward in the disk rather than 
inward, or a modified form the of the magnetic flares of 



Poutanen & Fabian (1999). Even if a scattering corona 
intercepts a significant fraction of the photons from the 
disk, negative phase lags could be observed so long as 1) 
the scattering time in the corona is much smaller than 
the intrinsic time lags (Miller 1995) and 2) the differ- 
ences between the energies of the seed photons and the 
observed signals are not too large (which should be the 
case in GRS 1915+105, since the apparent disk tempera- 
ture is > 0.8 keV; Nowak & Vaughan 1996). 

Moreover, the radius of the Compton corona can change 
without greatly affecting the slope and cut-off energy of the 
power law. The electron temperature (kT e ) of the corona 
is constrained to be nearly constant by the cut-off in the 
power law at 70-80 keV measured during X-ray faint ob- 
servations (left and center panels of Figure 3) . If the tem- 
perature of the corona is fixed, a steeper power law (higher 
T) implies a lower optical depth (see Sunyaev & Titarchuk 
1980). Since the power law is steeper during radio plateau 
observations (r ~ 2.5) than during the radio faint, X-ray 
faint observations (r ~ 2.1), we can suggest that the opti- 
cal depth of the corona decreases and its radius increases 
as the compact radio jet becomes more luminous. How- 
ever, the X-ray bright, radio faint observations would then 
have a small, optically thin corona, since the phase lags 
are negative and the power law is steep (r ~ 3.2) during 
these observations (Figures 3 and 4). The electron accel- 
eration mechanism could be different in the X-ray bright 
hard state, as is hinted at by the absence of an observable 
cutoff in the power law below 200 keV. 

As others have noted, estimates of the Lorentz factors 
for electrons emitting synchrotron emission (~ 400; Falcke 
& Biermann 1999) suggest energies that are much larger 
than that of the Comptonizing electrons. One may con- 
sider the acceleration mechanism for the jet as an unde- 
termined structure that draws material from the Compton 
corona near the inner disk, as in Fender et al. (1999a). In 
the following two sections, we therefore briefly speculate 
on the implications which the simplest geometries for the 
Compton corona described in the literature have on the 
jet production mechanism. 

4.2. A Spherical Corona within the Inner Accretion Disk 

Many authors have suggested that a spherical corona is 
formed within the inner radius of an accretion disk, either 
as an advection-dominated accretion flow (ADAF; Esin, 
McClintock, & Narayan 1997), or by a shock in the accre- 
tion flow (Chakrabarti & Titarchuk 1995). The compact 
jet responsible for the radio emission could form from a 
wind emanating from an ADAF (Blandford & Begelman 
1999) or a post-shock region (Chakrabarti 1999). The dis- 
tinctive feature of the spherical corona model is that the 
radius of the inner accretion disk is variable, which leads 
naturally to a scenario in which the X-ray timing and spec- 
tral properties, and the strength of the radio emission, are 
set by the inner radius of the accretion disk. 

The X-ray timing property that is perhaps most relevant 
to discuss in the spherical corona model is the 0.5-10 Hz 
QPO. The QPO may represent a characteristic time scale 
in the accretion flow such as an acoustic (e.g. Chen & Taam 
1994; Abramowicz, Chen, & Taam 1995), Keplerian (Ro- 
driguez et al. 2000), or free-fall (Molteni, Sponholz, & 
Chakrabarti 1996) time scale. Since the QPO peak is nar- 
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row {Q > 3), it must originate in a localized region, as all 
of the time scales increase as the radius increases. Among 
the possibilities for the origin of such a QPO are oscil- 
lations in a shock in the accretion flow (Chakrabarti & 
Manickam 2000) or a spiral density wave in the inner disk 
(Rodriguez et al. 2000). Since the frequency of the QPO 
decreases during episodes in which the inner radius of the 
disk is thought to increase (Belloni ct al. 1997; Markwardt 
et al. 1999), it seems plausible that the frequency of the 
0.5-10 Hz QPO in GRS 1915+105 tracks the inner radius 
of the accretion disk. Figure 7a then implies that the radio 
emission is strongest when the inner disk is farthest from 
the source. 

Nobili et al. (2000) have developed a model for the phase 
lags in GRS 1915+105 by assuming a spherical corona 
with two temperature regions: a warm (T ~ 1.5 keV) 
outer corona and a hot (T ~ 15 keV) inner corona. The 
corona is isotropically illuminated by soft seed photons. 
The sign and magnitude of the phase lags are correlated 
with the radius of the inner accretion disk, because as the 
disk moves inward it compresses the corona, increasing its 
optical depth. Large, positive phase lags occur when the 
inner disk radius is large, because the hot corona has an 
optical depth th ~ 1, while the warm corona is optically 
thin. When the inner radius is a factor of ~ 3 smaller, 
the optical depth of the hot corona becomes extremely 
high (th > 100), and the soft seed photons from the 
disk thcrmalize to the temperature of the optically thick 
hot corona. The thcrmalizcd photons are then Compton 
down-scattered in the cool corona, which produces soft 
lags. This model does not attempt to explain the phase 
lags as a function of Fourier frequency, but it does predict 
that if the QPO frequency tracks the radius of the inner 
disk, negative phase lags are expected when the QPO fre- 
quency is high — exactly as Reig et al (2000) observed in 
GRS 1915+105. However, a corona with an optical depth 
of t > 100 would produce a thermal spectrum, while the 
energy spectrum of GRS 1915+105 between 25-200 keV 
can not be described by an extremely optically thick spec- 
trum (e.g. a spherical Compton cloud with optical depth 
t > 5; Sunyaev & Titarchuk 1980) when the phase lags 
are negative. Further work needs to be done to reconcile 
the negative phase lags with the spectral shape under this 
type of model. 

The interpretation of the coherence function and the 
continuum power as a function of energy is less certain. 
When the QPO frequency is low, the coherence at low fre- 
quencies is nearest to (Figure 7c) and the ratio of the 
low frequency continuum power at high to low energies is 
near 1 (Figure 7d). A coherence less than 1 can occur if 
the variability is the superposition of separate linear sig- 
nals with different properties (Vaughan & Nowak 1997). 
When the inner disk is farthest from the source (low QPO 
frequency) there may be more area available for multiple 
inputs to contribute to the variability. 

The radius of the inner disk should also determine the 
energy spectrum under a spherical corona model. We have 
already discussed the model of Nobili et al. (2000). A dif- 
ferent set of assumptions are used by, for example Esin et 
al. (1997) and Chakrabarti & Titarchuk (1995), who take 
into account the thermal balance between the soft photons 
and the Comptonizing electrons. Soft spectra are formed 



when the inner radius of the disk is close to the last stable 
orbit, so that the disk produces copious amounts of seed 
photons which cool the corona. Hard states are formed 
when the inner radius of disk is far from the last stable 
orbit and most of the corona is photon-starved. However, 
if we consider the faint X-ray observations, the power law 
is steeper (r ~ 2.5) when the QPO frequency is lower and 
the disk is inferred to be farther from the black hole (ra- 
dio plateau observations; Figure 3 and 4) . This is opposite 
the sense expected, because the Compton corona should 
be cooled less efficiently (which implies a lower T) when 
the inner disk is large. If the QPO frequency tracks the 
transition radius between the disk and the corona, there 
must an additional factor beside the radius of the inner 
accretion disk which affects the observed X-ray properties 
of the corona in GRS 1915+105 under this model. 

4.3. A Planar Corona Sustained by Magnetic Flares 

In a planar corona model, one can assume that mag- 
netic flares powered by the differential rotation of the ac- 
cretion disk produce a population of relativistic electrons 
above the disk (di Matteo et al. 2000; compare Dove et 
al. 1997 for difficulties with static models). The compact 
jet could represent a collimated wind which forms accord- 
ing to the magnetic-centrifugal mechanism of Blandford 
& Payne (1982), and is strongest when the flares pro- 
vide electrons enough energy to escape the potential of 
the black hole (Romanova et al. 1998). In the magnetic 
flare model of di Matteo et al. (1999), the inner radius of 
the accretion disk in general does not vary. 

A natural starting point for discussing the planar corona 
models are the phase lags observed in GRS 1915+105. The 
time lags which would be implied by the scale height of 
the flares in the model of di Matteo et al. (1999) are 
too short (10 7 cm = 10~ 3 s) to explain the phase lags 
observed in black hole candidates. However, Poutanen & 
Fabian (1999) have developed a model for Cyg X— 1 in 
which the spectral evolution of the flares which accelerate 
Comptonizing electrons produce the observed power den- 
sity spectra and phase lags. Their model produces only 
hard (or zero) lags, because the electrons are assumed to 
be in thermal equilibrium during the flare, so the spectrum 
of the flares is set by the ratio of the heating in the flares 
to the amount of soft seed flux impinging on the flares. 
The flares become harder because the majority of the soft 
flux is produced by feedback in response to the flare, and 
the feedback is taken to decrease as a function of time as 
the flares move away from the disk. However, if the elec- 
trons accelerated in the magnetic flare have a non-thermal 
energy distribution which evolves as a function of time, it 
is conceivable that either hard or soft phase lags could be 
produced. The hard phase lags may represent the most 
efficient acceleration of the electrons, which could in turn 
explain why the radio plateau emission is associated with 
the largest hard phase lags. On the other hand, the nega- 
tive phase lags could represent spectral evolution in which 
the soft energy band peaks later, and relativistic electrons 
are accelerated less efficiently. 

In the case of a magnetically sustained planar corona 
with a fixed inner disk radius, there is some difficulty in ex- 
plaining what sets the frequency of the QPO and why this 
frequency varies. Since the power in the X-ray variability 
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in GRS 1915+105 drops sharply at frequencies higher than 
the 0.5-10 Hz QPO (Figure 4), we can speculate this QPO 
represents a resonance at the minimum time scale at which 
magnetic flares can form. The radio emission, then, would 
occur most favorably when only relatively slow flares form 
(low QPO frequency). 

The coherence function and the ratio of the high to low 
energy broad-band power would also be set by the evo- 
lution of the magnetic flares. Since the coherence at low 
frequencies is lowest when the QPO frequency is low and 
the phase lags are positive, the long, hard flares may un- 
dergo significant non-linear spectral evolution. The ra- 
tio of broadband power in the high and low energy bands 
could be set by the shape of individual flares in the respec- 
tive energy bands. Mathematically, this could be described 
in a manner similar to the shot noise models of Shibazaki 
et al. (1988). 

The model of di Matteo et al. (1999) predicts the en- 
ergy spectrum from magnetic flares, by assuming that a 
fixed fraction of the accretion power is released into the 
magnetic structures, and that the electrons which are ac- 
celerated by the flares are cooled inverse- Compton scat- 
tering of both cyclo-synchrotron and thermal (disk) seed 
photons. Soft spectral states are caused by flares which oc- 
cur close to the disk and are flooded by seed photons, while 
hard states are caused by photon-starved flares which oc- 
cur high above the disk. Although di Matteo et al. (1999) 
do not attempt to explain the phase lags observed in black 
hole candidates, we might expect hard phase lags to cor- 
respond with flares that produce the hardest (lowest T) 
power law, since the flares would be highest above the 
disk and cooled least efficiently. However, when the phase 
lags are positive in GRS 1915+105 we find a steeper power 
law (r ~ 2.5 in the plateau state) than when the phase 
lags are negative (r ~ 2.1 during radio faint, X-ray faint 
observations). Clearly, simple assumptions about the evo- 
lution of magnetic flares would have to be relaxed in or- 
der to explain the X-ray timing and spectral properties of 
GRS 1915+105. 



5. CONCLUSIONS 

We have found that although radio emission is always 
associated with the hard-steady state of GRS 1915+105 
(between 1996-1999), the 15.2 GHz radio flux density and 
the 2-200 keV X-ray flux are not correlated. Instead, we 
find that the hard states can be described using three ba- 
sic regimes of radio emission: a) radio plateau conditions 
in which bright, optically thick emission originates from a 
compact jet, b) radio steep conditions in which bright, op- 
tically thin emission originates from material most likely 
decoupled from the X-ray emitting region, and c) faint ra- 
dio conditions for which our current study is unable to 
constrain the spectrum of the radio emission. As the ra- 
dio flux increases, 1) the frequency of a ubiquitous 0.5- 
10 Hz QPO decreases, 2) the Fourier phase lags between 
hard (11.5-60 keV) and soft (2-4.3 keV) photons in the 
frequency range of 0.01-10 Hz change sign from negative 
to positive, 3) the coherence at low frequencies decreases, 
and 4) the relative amount of low frequency power in hard 
photons compared to soft photons decreases. 

We have attempted to understand these results quali- 
tatively using simple models collected from the literature, 
but we find that our attempts to understand the X-ray 
timing properties of the hard state of GRS 1915+105 in- 
evitably lead to contradictions with the observed energy 
spectrum. These observational results therefore provide 
impetus for developing more sophisticated models of the 
accretion flow and the compact jet in GRS 1915+105. 
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Fig. 1. — Monitoring data from X-ray and radio wavelengths. Top panel: Count rate as a function of time from the RXTE ASM, 1.5—12 
keV. Vertical bars across the top represent hard-steady pointed PCA/HEXTE observations. Dots represent the remainder of the pointed 
observations. Second panel: RXTE ASM hardness ratio HR2 (5-12 keV / 3-5 keV). Third panel: Flux as a function of time from the Ryle 
telescope at 15.2 GHz. Bottom panel: GBI radio spectral index (a = A log S v /A log v, where A log S v is the difference between the logarithm 
of the mean daily flux at v = 8.3 and 2.25 GHz; points are plotted only if the flux density at both frequencies was larger than 20 mjy). 
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Fig. 2. — The radio flux density at 15.2 GHz measured with the Ryle telescope as a function of the 2-200 keV X-ray flux measured with the 
PCA aboard RXTE. All hard-steady observations exhibit detectable radio emission, although the strength of the X-ray and radio emission 
are not correlated. The radio conditions of each observation are indicated as follows: triangles represent radio-plateau hard-steady conditions, 
circles radio-steep conditions, and x's radio-faint conditions. Three characteristic observations are indicated in the diagram: (1) a radio-faint, 
X-ray faint observation, (2) a radio plateau, X-ray faint observation, and (3) a radio-faint, X-ray bright observation. 
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Fig. 3. — The energy spectrum and the corresponding spectral fits, as described in the text, for three characteristic observations: a radio- 
faint and X-ray faint observation on MJD 50488, a radio-plateau observation on MJD 50737, and a radio-faint and X-ray bright observation 
on MJD 50708. The power law component dominates the spectrum at all energies, while the multi-temperature disk and the Gaussian iron 
line contribute a smaller fraction of the flux. 
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Fig. 4. — Timing properties of the three characteristic observations in Figure 3. Top panel: PDS in two energy bands, the black line 
indicates the 2-4.3 keV energy band, while the grey line indicates the 11.5-60 keV band. Middle panel: the phase lags between the 2-4.3 kcV 
and the 11.5-60 keV bands. Bottom panel: the coherence function between the 2-4.3 keV and the 11.5-60 keV bands. In all panels, we have 
indicated the frequencies of the QPOs in our fits to the 2.5-5 keV band with black short-dashed lines, and the break frequencies with black 
long-dashed lines. We have also indicated the break frequency in the 11.5-60 keV band with grey long-dashed lines (see below). 
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Fig. 5. — The relation of the timing properties to the X-ray energy spectrum for the three observations in Figure 3. Top panels: The RMS 
variability at all frequencies as a function of energy plotted with diamonds connected by the solid line. Middle panels: The mean phase lag 
between the low frequency break and the QPO frequency (region 2) as a function of energy plotted with diamonds connected by the solid 
line. Positive values indicate that the hard photons lag the soft photons. Bottom panels: The average coherence function measured below the 
low frequency break (region 1) plotted with diamonds connected by the solid line. All Panels: Fraction of the count rate contributed by the 
multi-temperature disk model (as derived from the spectral fits) as a function of energy, plotted as triangles connected by the dashed line, 
and labeled on the axes to the right. 



Radio and X-ray Emission from GRS 1915+105 



17 




Fig. 6. — The time evolution of the parameters of the PDS and CS for a stretch of faint radio emission (MJD 50400-50600, left panel) 
and bright radio emission (MJD 50725—50755, right panel), a) The ASM light curve, 1.5-12 keV. b) Ryle monitoring data (15.2 GHz), c) 
The radio spectral index (a) measured by with the GBI (see Figure 1). d) The 0.5-10 Hz QPO frequency (measured from the 2-4.3 keV 
band), c) Phase lags for the 2—4.3 keV relative to the 7.8—11.5 keV band (diamonds connected with the solid line) and the 11.5-60 keV band 
(triangles connected with the dashed line) measured between the low frequency break and the QPO. f) The average coherence function for 
the 2—4.3 keV relative to the 7.8—11.5 keV band (diamonds connected with the solid line) and the 11.5—60 keV band (triangles connected with 
the dashed line) measured below the low frequency break, g) The RMS power at frequencies below the QPO measured in the 2-4.3 keV band 
(diamonds connected with the solid line) and in the 11.5-60 keV (triangles connected with the dashed line). 
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Fig. 7. — Parameters of the PDS and CS as a function of the radio flux (15.2 GHz). Triangles represent radio-plateau hard-steady 
conditions, circles radio-steep conditions, and x's radio-faint conditions, a) The QPO frequency, b) Phase lags for 11.5-60 keV photons 
relative to 2-4.3 keV photons, averaged over intermediate frequencies (between the low frequency break and the QPO). c) The coherence 
function for 2-4.3 keV photons relative to 11.5-60 keV photons, averaged over low frequencies (below the low frequency break), d) The ratio 
of the integrated broad band powers below the QPO frequency in the high (11.5-60 keV) to low (2-4.3 keV) energy bands. 
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Fig. 8. — Correlations between the QPO frequency (top panels) or the average phase lag at intermediate frequencies (11.5-60 keV relative 
to 2-4.3 keV; bottom panels) with fluxes derived from spectral fits to the data. Left panels: total flux. Center panels: thermal flux, from the 
multi-temperature disk model. Right panels: power law flux, integrated from 2-200 keV. The symbols are the same as Figure 7. 
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Table 1 

Hard-Steady Pointed RXTE Observations of GRS 1915+105 



ObsID 




Date 


X-ray 


HR b 


$2 25 


$ 8 3 


5*15 2 






(UT) 


Flux a 




(mjy) 


(mjy) 


fmJv) 


10408-01-22-00 


1996 


Jul 11 02:09 


4.15 


1) 


0.063 






54(2) 


10408-01-22-01 


1996 


Jul 11 05:16 


3.93 


1) 


0.070 






54(2) 


10408-01-22-02 


1996 


Jul 11 08:38 


3.78 


1) 


0.073 






54(2) 


10408-01-23-00 


1996 


Jul 14 11:45 


4.17 


1) 


0.063 






39.3(6) 


10408-01-24-00 


1996 
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a Between 2-200 keV, in units of 10~ 8 erg cm" 2 s" 1 . 
b Ratio of PCA counts from 12-60 : 2-12 keV. 

Observations which were not analyzed either because they were not taken with the 
necessary data modes or were too short for meaningful timing analyses. 
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Table 2 

Radio Observations of GRS 1915+105 while Faint 



Date Time Instrument Frequency S v Beam 
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35 - 26 


0.63 


x 0.36 


1997 


Jan 12 


21:50 


- 22:30 


VLA 


8.4 


30 - 19 


0.41 


x 0.29 


1 QQ7 


Tan 1 9 

o all LZi 


12:45 


- 20:57 


GBI 


2.25 


44(4) 




11 


1997 


Jan 12 


12:45 


- 20:57 


GBI 


8.3 


21(6) 




3 


1997 


Jan 13 


16:50 


- 18:35 


VLA 


8.4 


4.5(1) 


0.27 


x 0.24 


1997 


Jan 13 


16:50 


- 18:35 


GBI 


2.25 


14(4) 




11 


1997 


Jan 13 


16:50 


- 18:35 


GBI 


8.3 


4(6) 




3 


1997 


Apr 19 


10:40 


- 15:20 


VLA 


1.4 


3.9(1) 


4.5 


x 3.5 


1997 


Apr 19 


10:40 


- 15:20 


GBI 


2.25 


6(4) 




11 


1997 


Apr 19 


10:40 


- 15:20 


GBI 


8.3 


6(6) 




3 


1998 


Sep 14 


04:20 


- 05:30 


VLA 


8.4 


8.2(1) 


0.74 


x 0.70 


1998 


Sep 13 


23:07 


- 04:26 


GBI 


2.25 


18(4) 




11 


1998 


Sep 13 


23:07 


- 04:26 


GBI 


8.3 


11(6) 




3 


1998 


Sep 13 


21:39 


- 22:39 


Ryle 


15.2 


6.0(6) 




30 


1998 


Sep 29 


04:00 


- 06:30 


VLA 


5.0 


12.2(1) 


2.26 


i x 1.2 


1998 


Sep 29 


04:00 


- 06:30 


VLA 


15.2 


13.9(1) 


0.40 


x 0.38 


1998 


Sep 28 


21:36 


- 03:36 


GBI 


2.25 


20(4) 




11 


1998 


Sep 28 


21:36 


- 03:36 


GBI 


8.3 


12(6) 




3 


1998 


Sep 28 


19:07 


- 21:30 


Ryle 


15.2 


10.7(4) 




30 


1998 


Oct 08 


03:30 


- 05:30 


VLA 


5.0 


31.7(1) 


1.6 


x 1.3 


1998 


Oct 08 


03:30 


- 05:30 


VLA 


15.2 


18.5(1) 


0.42 


x 0.41 


1998 


Oct 07 


21:22 


- 03:07 


GBI 


2.25 


60(4) 




11 


1998 


Oct 07 


21:22 


- 03:07 


GBI 


8.3 


23(6) 




3 


1998 


Oct 07 


19:05 


- 21:19 


Ryle 


15.2 


21.3(4) 




30 



